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Abstract 

The interaction of orbital electrons with the charge and magnetic moment of the nucleus polarizes it, 
and the detailed description requires a careful treatment of the nuclear vector polarizability. We present 
a complete and closed form expression for the resulting contribution to hyperfine splitting in light atomic 
systems. 
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Atomic energy levels are split due to the interaction between the nuclear and electron magnetic 
moments, /2 and jle respectively. This splitting, called the hyperfine splitting (hfs), for a non- 
relativistic hydrogen-like system in an S'-state is given by the Fermi contact interaction (in units 

h = c = 1,6"^ = An a) 
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\lI-fIe6%rm = ^^gS.a 
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where 



Ze 
2M 
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g = {1 + k)/S, and Z e and M are the charge and mass of the nucleus. For electrons the (^-factor is 
close to 2, with a small anomaly k ^ a/(2 tt), which is neglected in Eq. ([U), and for the nucleus we 
assume an arbitrary (^-factor. To compare with high accuracy measurements of hyperfine splitting, 
such as in the deuterium atom = E/h) 



z/exp = 327384.352 522 2(17) kHz, 



(3) 



many higher order relativistic and quantum electrodynamics (QED) corrections have been calcu- 
lated with high precision |2]. However, the theoretical accuracy is limited by nuclear effects. For 
example, in deuterium yi |41] 
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138 ppm, 
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where z^qed is the QED prediction assuming a point nucleus. This difference is many orders of 
magnitude larger than the precision of z/cxp and i'qed, and it is attributable to the nuclear structure 
correction. We note that the leading order 0{Z a) relativistic correction vanishes for a point-like 
and infinitely heavy nucleus. The subject of this work is the detailed study of the effect of both 
finite nuclear mass and nuclear structure, which make this 0{Z a) correction nonvanishing. This 
correction can be represented by the two-photon exchange forward scattering amplitude, which in 
the temporal gauge A° = takes the form (cu = 



6E, 



hfs 
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where 
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p is the momentum at rest, T^^ is the corresponding virtual Compton scattering amplitude off the 
nucleus, and a subtraction of the linear divergence at small k in Eq. ([5]), which is related to the 
leading Fermi interaction, is assumed implicitly. 

One splits the nuclear structure correction 5i?hfs into three parts: 

SEi^ow is the leading correction to hfs of order Z amrN, where is the size of the nucleus. This 
correction has been first derived by Low in Js], and has been recently reanalyzed and calculated 
for such nuclei as D, T, and '^He by Friar and Payne in Ref. [3] . In this work we present a complete 
derivation of the Low correction, as well as the higher order Z a m/nip correction which comes 
from nuclear excitations and recoil, which we denote by SEp^h The second term in Eq. (|7]), 
5-E'zomach is the SO Called Zemach correction from individual nucleons. For the hydrogen atom it is 



^-Ezemach 2am f (Pk 



Eus vr2 7 k^ 



GEi-k^)GMi-k^_^ 



—2amrz, (8) 



1 + K 

where Ge and Gm are the electric and magnetic formf actors of the proton. It is convenient to 
rewrite this correction in terms of the Zemach radius rz 

'^z = y (frpEir) pM{r) r, (9) 

with pe and pM being the Fourier transforms of G^; and Gm/{^ + f^)- The results of Ref. yfl 
for the proton and neutron Zemach radius are 1.086(12) fm and —0.042 fm respectively. For an 
arbitrary nucleus this correction is a coherent sum of Zemach corrections from all nucleons, 

(^-E^Zemach = ^ ^^^-^ • (-2am) {y^^gaSaTaz)- (10) 
^ a 

One notes that it would be not accurate to only use the Zemach formula with the nuclear elastic 
formfactors. Instead one should assume a point nucleus for the QED part of hfs^and at the first 
approximation calculate ^-Elow and (J-Ezemach as was done for example in Ref. [3]. It is still an 
open issue as to the accuracy of the elastic formfactor treatment of hfs in heavy atoms or ions, and 
its relation to SEi^o^, but this problem is not studied here. 

We derive below a complete closed form expression for 5Epoi which can be used to improve 
theoretical predictions for hyperfine splitting in light atoms or ions. The main idea behind this 
derivation is the existence of an expansion parameter in the effective nuclear Hamiltonian, namely 
the ratio of the characteristic momentum Q of a nucleon to its mass nip, which is about 0.10 — 0. 15 



in typical nuclei IJ]. ^-E'poi accounts for all the nuclear structure corrections, which are Q/mp 
smaller than the leading ^-Elow contribution. To carry this out we split the integral in Eq. ([5]) into 
two parts. In the low energy part, where k is of order of the binding energy per nucleon, both 
the elastic contribution and nuclear excitations play a role. In the high energy part, where k ^ 
binding energy, nucleons are seen by the electron as free particles, and their binding energy can 
simply be neglected. This derivation will be similar to the one presented by Khriplovich et al. in 
Refs. Da lZD for the particular case of the deuterium atom, but differs in many details and in the 
resulting formula. In particular we show that the leading logarithmic contribution vanishes, while 
the result of Ref. [|7|] can not be rewritten to such form. 

We start the derivation by noting that according to Eq. ([5]), only the antisymmetric part of the 
Compton amplitude T*^ contributes to hfs. If the scattered photon is on mass shell, uj = 1^1 > then 
this amplitude T*^ = iuj"^ a*-' can be expressed in terms of the vector polarizability d"^ [|8|]. In the 
general case uj ^ |A;|, and for convenience we will identify the antisymmetric part of T*-' with the 
vector polarizability 

= {^'^ -T^')I1. (11) 
Then using Eq. (|6l), the nuclear structure correction to hfs is 



hfs 
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[J^ e^^^ + k^ e'*^' - U k^ e^^^) a- 
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(2 Tif UJ (cj2 _ A;2) (cu2 - 2 m - fc2) (cj2 + 2 m - fc2) 

where the Feynman integration contour is assumed and the apparent 1/uj singularity cancels out 
with the numerator. For small photon momenta the elastic part of the vector polarizability 
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can be obtained from the following Hamiltonian, which describes the interaction of a particle 
having charge Z e, mass M, and spin S with the electromagnetic field 

n2 



H 



Ze 



T79S-B 



Ze 



{g-l)S -[E xll-llx E]. 



(14) 
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where n is defined in Eq. (fTSl ). Eq. (fT3] ) is in fact a low energy virtual Compton scattering 
amplitude off a point nucleus with arbitrary spin S. It agrees with results obtained for the first 
time by Khriplovich et al in Ref. [6]. To obtain the corresponding contribution to the hyperfine 
splitting (^o-Ehfsj we assume \k\ < A, with the cut-off A being larger than binding energy per 
nucleon, and smaller than the inverse of the nuclear size. We will show later that the dependence 
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ln( (/-4(7-12) + ^((?2 + l24(7 + 4) 
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on A cancels out between all the corrections, so there is no need to specify its precise value. After 
subtraction of the leading nonrelativistic part ~ A;^, which is a Coulomb iteration of the Fermi 
contact interaction, and after the to and k integrations of ^qE^^s are performed, it takes the form 

The important point of this calculation is that after expansion in the small parameter m/A or 
AE/A, the limit A ^ is performed consistently in all the parts. 

Apart from the elastic contribution, various nuclear excitations play a role, even for small values 
of the photon momentum k. Their calculation is more complicated, as the nuclear Hamiltonian 
is not well understood. We will use the interaction Hamiltonian with the electromagnetic field 
which was obtained by us in [9]. The main assumption is that the characteristic wavelength of 
the electromagnetic field is much larger than the nuclear size. Although it was derived for a 
system consisting of electromagnetically interacting particles, we assume that the obtained form 
should be valid also for nucleons: this follows from the gauge and Lorentz symmetries of the 
nucleon-nucleon interactions. However, this assumption can be explicitly verified only by using a 
systematic method to derive the nuclear Hamiltonian such as chiral perturbation theory (xPT), but 
this will not be done here. This Hamiltonian, dropping terms that do not contribute to the vector 
polarizability, is \Q\ 
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where the index a goes over protons and neutrons, M = J2a ""^a and Z e = e^. For protons 
Ca = e, ga = gp, while for neutrons eaga = e 9n, but alone 0. The electromagnetic field and 
its derivatives are assumed in the above equation, to be at point R. Center of mass coordinates are 
defined by 



R 

n 



— m„ fa , 

eaA{R) = P-ZeA{R) . 



a 



(17) 
(18) 



The relative coordinates, which are defined by 



Xa = fa-R, (19) 

ma 
M 



qa=Pa-'^P, (20) 



K,Qi] = ^S^5^,-^]. (21) 



obey nonstandard commutation relations 

1 = i i)'^ ( ^.u - 

M 

The total spin S is 

S = J2^a + l'a, (22) 
a 

where the orbital angular momentum is 

la = XaX Qa, (23) 

and the ^f-factor is defined by 

a 

i^iN is the internal Hamiltonian of the nucleus, the exact form of which is not well known. For 
further calculations we will assume that the following commutator 

[Hi^,x\] = -z^, (25) 
rrir, 



"a 



holds, at least to a good approximation. 

Since the nuclear excitation energies are larger than the electron mass m (with a few excep- 
tions), we neglect m in Eq. (fT2l) and obtain 

There are various contributions to a'''' which follow from the Hamiltonian in Eq. (fT6l ). and we 
calculate them in order. The first contribution is due to the electric dipole coupling 

6H = -Y,eaXa-E = -D-E. (27) 

O 

This contribution has already been considered by Friar and Payne in Ref. mOD for the particular 
case of the deuterium atom. For the general nucleus the effect on the vector polarizability is 
analogous 



and the corresponding contribution to the hyperfine splitting is 

3i 



16 TT^ 



(29) 



The constant m (the electron mass) in the argument of the logarithm is not relevant here, as it does 
not affect the above matrix element, since D' commutes with . In the nonrelativistic approxi- 
mation ^i-Ehfs vanishes. Only the presence of relatively small spin-orbit terms in i/jN causes it to 
be nonvanishing, therefore this correction (JiE'hfs is not expected to be the most significant one. 
The second contribution is due to magnetic dipole coupling in Eq. (fT6l) 



6H 



2mr, 



The corresponding vector polarizability is 

/ r/^ rik 
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and the contribution to the hyperfine splitting is 

1 



[ill - ill)) X k\' 
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where 



^^,:|ff.fll = E — 

M ^-^ rria 

a 

The third contribution is due to the electric quadrupole coupling 



sh = -Iy1 « - '^^V3) e:^ ^ -1 D^^ Ey . 



(30) 
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The corresponding vector polarizability is 
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and the contribution to the hyperfine splitting is 
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where 



M 



^ i [H, D'^] = J2^[xiqi + qi x{ - (f „ + q, f„) 5^73] . 



(37) 



The contribution coming from the magnetic dipole on one side and the electric quadrupole on 
the other side is being neglected. This is because the resulting matrix element 

{Q'^' ln[(i7iN - E)lk\ / - h.c ), (38) 

involves operators Q^^ and ^i^ which commute when the implicit sum over k is assumed. Therefore 
this matrix element does not depend on A, is small, and will consequently be neglected. This 
argument will be used several times in neglecting or simplifying expressions in the following. 
The next contribution comes from relativistic corrections to the electric dipole coupling. 

The vector polarizability correction is 



{9a - 1) 
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and the contribution to the hyperfine splitting is 
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We have used in the above an approximate relation 

{Q" ln[2 (//in - E)/m] < - h.c. ) ^ (Q'^ ln[2 (Hj^ - E)/m] x^ < - h.c. )/3, 

which comes from the fact that the commutator [Q^ , x\ — x^ 8^/3] vanishes. 
In the fifth contribution one vertex remains -DE but the other one is 



(41) 



(42) 
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The vector polarizability is 

^ 2ma 
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Xak (Sq X A)) 
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Xak {Sa X k)'' — {k l,ljJ 
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and the corresponding correction to the hyperfine splitting is 

^'■^'(0) 17 E ('a^i ln[2 (i/iN - E)/m] - h.c.) 



I e 
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The sixth contribution is due to the following modification of one of the vertices 

' ' in m ' 



(46) 



It is very similar to the previous one, and can be obtained by replacing spin with the orbital angular 
momentum. The vector polarizability is 
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and the corresponding correction to the hyperfine splitting is 
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The seventh contribution is due to a relativistic, spin independent correction to the electric 
dipole operator, 

^H = Y,^ X - X /;) ■ dtE. (49) 

a 

The vector polarizability can simply be obtained from Eq. (|40|) 



-2 UJ 
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+ (/(! X Xfl Xq, X /(j) 



E-Hw + u 
and the contribution to the hyperfine splitting is 
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We have used in the above the approximate relation 

<E « + <) ln[2 (//in - E)/m] - h.c.) 

a 

^ (E + ^'«) (//in - /?)/m] - h.c.>/4, 



(52) 



which comes from the fact that the commutator [Q^ , (x^ + l{ x^) — (x{ + l\ x{)/4] vanishes. 

The remaining two low energy contributions involve the momentum of the nucleus P. The 
following correction 



5H 



D 



2M 



-(P X B - B X P). 



(53) 



replaces the electric dipole coupling —D E. The resulting vector polarizability is 
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— {k ^ /, a; — >• —a;), 
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and the contribution to hyperfine splitting is 

3e2 



6^E, 



^a iV^'(0)^(^xln[2(//iN-^)/m]g + ln[m/(2A)]Dx g). (55) 



The kinetic energy of the nucleus 



5H 



P^ 
2M'' 



(56) 



modifies the vector polarizability in Eq. (|28l) by adding A;^/(2M) to //in which results in the 
following correction to the hyperfine splitting: 



SoE, 



9-C/hfs 



^ no) ^{Dx ln[2 (//in - E)/m] Q + 



, , m \ 29 
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The last contribution comes from photon momenta \k\ > A. Since A is much larger than the 
binding energy per nucleon one can completely neglect the nucleon-nucleon interaction. There- 
fore, the electron effectively sees free protons and neutrons. In the case, when a photon is emitted 
and absorbed by the same nucleon, we replace 

r-A 

dk{...) (58) 

and separately consider both terms. The second term gives a contribution similar to SoE^is, but 
with the opposite sign: 



/ dk{.. 
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(59) 
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while the first term of Eq. (1581) consists of a Zemach correction which is already included in Eq. 
([7]) as 5-Ezemach, and a small polarizability with recoil corrections from individual nucleons, which 
we account for in 5g 

SiAs = y,7r^ ^\'^)asa5ga. (60) 

a 

For the proton the sum of the recoil and polarizability corrections, 5.84 ppm and 1.30 ppm respec- 



tively, amounts to llllh 



^ = 7. 14 ppm, (61) 

while for the neutron it has not yet been obtained, although we expect it to be much smaller. 

We now consider the case where the photon is emitted and absorbed by different nucleons, 
denoting the resulting correction as (5-Elow The Compton amplitude T^^ obeys T'^" = 0. From 
this one obtains /c* T*'^ = tuT'^'', and uses this to rewrite the formulae (|26l) in the form 

SE,,. = ^.^(0) / If / (||; + (T»' - r»)] . (62) 

Since the dominating contribution in the above integral comes from uj,k of order of the inverse 
of the nuclear size, which is much smaller than the nucleon mass, we may use the nonrelativistic 
approximation for T'^'^, 



r.-1-V. J \ 



L — /iiN — uj + te 



+r,{r') e-^^^^ rr^^ ^ - J'air) e^''). (63) 
E - Hm + UJ + te / 

The contribution coming from T'^' in Eq. (|62|) is of nominal order Z am /nip, therefore after 
carrying out the cu integration one can neglect //in — E in comparison to k. It is then proportional 
to the commutator [J^ (r), Jl,{r)], which vanishes for different nucleons a ^ b. The second term 
in Eq. (|62|) involving is of nominal order Z amr^, therefore one shall keep also the second 
term in the expansion: 



1 1 Hw-E 



k + Hm-E k P 
However, this second term leads to the correction proportional to 



(64) 



(J°(r) e*^'^^(iJiN - E) Jlir') e"^^^^" + J^(r') e"^^^^" {H,^ - E) J'^{r) e'^^') 
([J°(r) e^^^n [i/iN - E, Jlir') e"^^^^"]]), (65) 
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which also vanishes for a^h. Therefore SEi^ow becomes 



<5ELow = e2V''(0) 



A (2 7r)3 fc6 



axkJ2 dhdh'{jl{r) e^'=(-^--") /,(/)>. (66) 



The linear divergence at small A is eliminated by subtraction of l/k'^ term from the integrand, 
which is related to the leading Fermi interaction. We should subtract it at the beginning, but for 
convenience we do it now. Assuming the nonrelativistic approximation for the interaction of the 
nucleon with the electromagnetic field, 

{p-eAf 
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2 m 



+ eA° - jl-B, 



(67) 



one obtains the "Low" nuclear structure correction to the hyperfine splitting [|50 
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A more accurate treatment of that includes the finite size of the nucleons is presented in 

n 

Ref. [3], where the authors emphasize the importance of additional meson exchange currents. 

We now present the complete nuclear polarizability correction to the hyperfine splitting 5-Epoi, 
and point out that it does not depend on the artificial parameter A 
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where Ln = In [2 {Hj^ ~ E) / m], m is the electron mass, Sg is defined in Eq. ( |60l ). -D* in Eq. (|27l ), 
/i* in Eq. (l30l) . in Eq. (l33l) . Q*-' in Eq. (l37l) . Moreover, our definition of the nuclear ^f-factor 
Eq. ^ differs from the standard one by the use of the actual nuclear charge and mass, as opposed 
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to the unit charge and the proton mass. The above expression is quite complicated, and its use 
may be limited by the lack of an accurate nuclear wave function. However, it takes a particularly 
simple form for the deuteron. 
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12 mm, 



(70) 



where R 



Tn and Ht and Hs are nonrelativistic proton-neutron Hamiltonians in the triplet 



and singlet states respectively. Further on, we used an approximate deuteron mass = 2 mp and 
neglect D-wave mixing of the ground state. 

The nuclear structure correction to hfs in deuterium, including some polarizability effects, has 
been obtained by Khriplovich et al. in Refs. |^ 7]. Their result of 153 ppm seems to explain 
well the difference between the experimental value and theoretical predictions with assuming a 
point nucleus, which is 138 ppm, see Eq. ([4]). We point out, however, that their calculations 
were performed only in logarithmic accuracy, and in our result the In A, the logarithm of an arbi- 
trary chosen cut-off cancels out. Friar and Payne have recently calculated in Ref. [|3|] SEi^^w and 
(^-Ezomach- Their result for deuterium of 141 ppm accounts for most of the difference in Eq. (31), 
but they have not included 5Epoi, which they claim to be negligibly small, in disagreement with 
our findings. It would be now interesting to calculate the nuclear polarizability correction to hfs 
according to Eq. (fTOl) and to verify it against precise experimental value for hfs in deuterium yj] in 
Eq. (|3]), as well as in other light atoms. As it is highly difficult to go beyond 5Epoi, the value of this 
correction sets the ultimate limit of the accuracy theoretical predictions for light atomic systems, 



which we expect to be a few ppm, as they are for hydrogen hfs [lllll . 



The nuclear polarizability correction 5Epoi is of order 0{Z am/nip), thus, in general, it is 
smaller than SE-^ow or SEzemach- However, in some cases we expect this correction to be more 
significant. When a closely lying state of the opposite parity is present, such as in ^^Be nu- 
cleus, with a large B(E1) ^ 0.1 e^fm^ line-strength [|l2n . SiE^fs in Eq. (|29l ) can be as large 



as B{E1) miripEp 10 ^ Ep, which should be possible to verify experimentally. In fact, the 
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measurement of the hyperfine splitting in ^^Be"*" has already been proposed in Ref. lll3ll for the 
determination of the neutron halo. We do not think that the contribution to hfs from the neutron 
halo can easily be identified, but this measurement can verify the significance of the nuclear vector 
polarizability in the hyperfine splitting of atomic systems. 
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